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I. INTRODUCTION 



o 

^jH" Rotating black holes in higher dimensions were first discussed in the seminal paper of Myers and Perry [1]. One 
^ \ surprising feature of the Myers-Perry (MP) solution is that in general j spin parameters are required. The first 
asymptotically non-flat five-dimensional MP metric was given in [2], where subsequent generalizations to arbitrary 
dimension, with multiple rotation parameters was done by Gibbons, Lii, Pope and Page (GLPP) [3, 4], which gave 
_ -. • a formal proof of the solution in [2]. The most general Kerr-(A)dS-NUT metric was found by Chen, Lii and Pope 

(CLP) [5]. These two papers represent a fundamental step to developing separable equations of motion. 
\q , Deriving wave equations for higher dimensional rotating black hole spacetimes [6] relies crucially on the method 
£\| ' of separation of variables from a Hamilton- Jacobi equation, see also [7]. Using the CLP metric, Frolov, Krtous and 
Kubiznak [8] were able to separate the geodesic equation and find the Klein-Gordon equation in the most general 
setting (this was generalized to charged cases in reference [9] for zero NUT parameter). The actual separation is due 
\Q \ to the presence of hidden symmetries in the form of Killing tensors [10-12], where a whole tower of Killing tensors and 
f^*) . symmetry operators [13] can be constructed with the help Killing- Yano and conformal Killing- Yano tensors. They 
guarantee the separability of the geodesic equation, the Klein-Gordon equation and also the Dirac equation [14]. 
Unfortunately the separation of the graviton does not appear to be possible based on these hidden symmetries [11]. 

In the literature research has been largely directed toward solutions with only one rotation parameter, the so 
called simply-rotating case, e.g., see [15]. However, recently, solutions with two rotation parameters in D > 6 have 
been investigated [16, 17]. In five dimensions, some work has been performed in the asymptotically flat case [6, 18], 
and in this paper we extend these results to non-asymptotically flat cases. These are relevant to the AdS/CFT 
correspondence, where black hole solutions in five-dimensional minimal gauged supcrgravity models require rotating 
solutions to avoid closed timelike curves (CTCs) [9, 19, 20]. The most general SUGRA solution in five-dimensions for 
arbitrary rotation parameters (ai, 02) was found recently in the work of Chong and collaborators [19]. The important 
point about the spheroidal harmonics in these cases is that they do not depend on the charge and hence will be the 
same as that for Kerr-(A)dS. 1 

This article is a follow up to [17] in which the Klein Gordon equation for doubly rotating black holes in D > 6 
dimensions was investigated. Superficially, D = 5 would not appear to be significantly different from D > 6; however, 
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1 Based on symmetry arguments in the metric we expect that the spheroidal harmonics are also independent of NUT charge. Hence our 
scalar spheroidal harmonics should actually be valid for charged Kerr-(A)dS-NUT black holes. 
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it turns out that in five dimensions only a single angular equation occurs, thereby making the analysis different in 
these two cases. 

In this work we will, after separating the Klein-Gordon equation on the GLPP background (with two rotations in 
five dimensions) investigate the (scalar) spheroidal harmonics, where in particular we show in detail how to apply 
perturbation theory to obtain the angular eigenvalues (separation constant) of the spheroidal harmonics. This is in 
the spirit of Fackerell & Crossman, and Seidel [21, 22], which used the properties of Jacobi functions. 

There are various reasons for studying this case. One is that in five dimensions the rotation parameter is bounded 
to something like a\UJ < 1.5 (for a 2 = 0) in units of black hole mass M = 1 and hence a perturbative expansion of 
the eigenvalue is well suited to speed up numerics for QNMs or greybody factors. Furthermore, as we mentioned, 
because there is only one spheroidal equation for a D = 5 doubly rotating black hole, the perturbative expansion is 
much simpler. 

The structure of the paper is as follows: In the next section (Section II), we discuss the general metric for Kerr- 
(A)dS black holes with two rotations in five-dimensions. The corresponding Klein-Gordon equation is separated into 
a radial and angular equation in Section III. In section IV we explain the perturbative method. Then in Section V (as 
a sanity check) we compare our analytics with the continued fraction method. Conclusions are then given in Section 
VI. 



II. CLP METRIC IN FIVE DIMENSIONS 



Many of the steps for the separation follow closely to that for D > 6 doubly rotating black holes [17], and hence, 
we shall only briefly outline the steps. As we are interested in the five dimensional case, we take the CLP metric to 
be 5]: 
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Here we shall only consider the Kerr-(A)dS case, so all NUT charges have been set to zero, where g 2 is related to the 
five dimensional cosmological constant with R^ u = —2>g 2 g^ u . Moreover, this metric is already in the Boyer-Lindquist 
form because there is no cross term for dr. The ingenuity of the CLP metric is the introduction of the coordinate y, 
which is related to the original direction cosines coordinates. In five dimensions these are fii = sin 8 and fi 2 = cos 9, 
with 
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y 2 = a\ cos 2 6 + a\ sin 2 9 . (2.13) 
This metric in equation (2.1) can be further simplified via the relations [5]: 
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In the above the linear relations are 
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The metric in equation (2.1) becomes 



rf.s 2 = ^dy 2 - ^dr 2 + ±1 [# - r 2 dVi] 2 + [#o + y 2 ^i] * 

Ax A2 Ui (72 

+ [^o + (2/ 2 - r 2 )d^ - y 2 r 2 d^ 2 ] 2 ■ (2.23) 
Note that, comparing with the CLP metric, see equation (22) in reference [5], we have 

4 0) = i , A^ = -r 2 , 4 0) = i , A 1] = y 2 , 

A<®=1 , A^=-r 2 +y 2 , A™ = -y 2 r 2 . (2.24) 



III. SEPARATION OF THE KLEIN-GORDON EQUATION 

Using the results of Frolov, Krtous, and Kubiznak [8] and considering their equation (4.2), the Klein-Gordon field 
in five dimensions has the following ansatz: 

$ = i? 1 ( 2 /)i? 2 (r)e l *°^ +l * 1 ' ,/ ' 1+4 * 2 ^ . (3.1) 

The separated equations for Ri and R 2 then follow: 
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where bi is the separation constant. 

By using the relationship between t, </>i, 4> 2 an d ipo, ipi, ip 2 hi Eqs. (2.20-2.22), we can obtain a relationship between 
the eigenvalues: 
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with < 9 < ir/2 or a 2 < y < a±, where b\ is a constant of separation. In the next section we shall effectively work 
with the variable y. 

To be able to compare our spheroidal harmonic with other work in the literature, it is also interesting to write an 
expression for the angular wave equation in terms of the latitude variable 9. From the angular equation in (3.2) above 
the derivative terms, with y 2 — a\ cos 2 9 + a\ sin 2 9, become 
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The non-derivative term can be simplified to: 



1 



1 



= -bi+LO 2 



V 



-bi 



2 2 9 2 



1 (l-ff 2 fl 2 )(l-ff 2 a 2 ) 
.9 2 9 2 (1- 9 2 (al cos 2 9 + a* sin 2 6)) 



+ w \—2{m\a\ + m 2 a 2 )] 



2/ i \2 '"l / 1 2 2\ i "'2 r-i 2 2\ 

g [nnai + m 2 a 2 ) + , 9 (1 - g %) H 5— (1 - g a 2 ) 



cos^ 



(3.8) 



Finally, the angular equation in terms of 9 becomes 
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This agrees with the angular equation in the general charged five-dimensional Kcrr-(A)dS case with two rotations 
that was found in [9], also see reference [6] for the flat g = case. 

For completeness, we also present the separation of the radial part, equation (3.3). We have 
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for the derivative terms, where we defined 

A = (1 + ,gV)(r 2 + a 2 )(r 2 + a 2 ) - 2Mr 2 . (3.11) 
In a similar way to the angular equation, we can simplify the non-derivative term in the radial equation as 
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This radial equation could be used to investigate QNMs, for example, on a five-dimensional Kcrr-(A)dS background. 



IV. EIGENVALUE EXPANSION FOR D = 5 



In Section III we have given a general analysis of the massless scalar equation in a Kerr-(A)dS spacetime with two 
rotations and D = 5 for general g. As we saw, the angular equation in this case has two rotation parameters, but 
only one angular equation, because there is only one latitude variable (cf. equation (2.13)). This simplifies things 
a great deal and hence we shall now work out the perturbative expansion of the angular eigenvalue using techniques 
similar to the works of [21, 22] (also see [23]). 

To develop a perturbative expansion for the eigenvalue b\ , it is convenient to make the change of variables 
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with —1 < x < 1. Then in terms of x, the angular equation becomes: 
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(4.2) 



where the constant B = b\ + 2uj{m\ai + m 2 a 2 ) — g 2 {rn\ai + m 2 a 2 ) 2 . Note that for simplicity we have dropped the 
subscript on R. 

Now comes the perturbative part of the method, which is essentially identical to time- independent perturbation the- 
ory in quantum mechanics [23]. Expanding in powers of a\ and a 2 , which we assume to be small, we can schematically 
write 
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The regularity of the form of the operators after O4 makes it possible to evaluate the eigenvalue iteratively as we shall 
do in the following. 

We start with the zeroth order in equation (4.3), 
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Before moving on we will also need to consider how the normalization condition above affects the higher order terms: 
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and using the result in equation (4.27), we have 

1 



(0 4 )m+4 = - — (a 2 -a 2 ) 2 g 2 uj 2 X lA+2 X l+2Mil (4.29) 
lb 

(0 4 )m+2 = -^(a\ - a 2 2 fg 2 uj 2 {Xu +2 X l+2 j +2 + X ul X ul+2 ) , (4.30) 

(04)i,i = ^(aj - a 2 ) 2 g 2 uj 2 {l - Xu +2 X l+2 + - X Ll X u - X M _ 2 X ( _ 2ii ) , (4.31) 

(0 4 )u-2 = -tt(o? - a 2 ) 2 g 2 u 2 {X U - 2 Xi- 2 ,i- 2 + X M X M _ 2 ) , (4.32) 
lb 

(0 4 )m-4 = - — (a?-ai)Vw 2 X M _ 2 Xi_ 2)/ _ 4 . (4.33) 



Putting these into equation (4.25), we have 

/mi 7712 

4 [(0 4 )m + (rf2)i ) / +2 (0 2 ) i+2 ,/ + (d 2 ) M - 2 (0 2 )/_ 2 ,/] 



1 2 2 2 f 8<? 2 w 2 [Z 4 + 4/ 3 + 2Z 2 (m 2 + m 2 ) + 4Z(m 2 + m 2 - 2) - 3(m 2 - m 2 ) 2 ] 
— (a 4 — a 2 ) 



64 v 1 " \ (/-l)J(f + 2)(Z + 3) 

[c 2 + g 2 l(l + 4)] 2 [(/ + 2) 2 - (\ mi \ + |m 2 |) 2 ][(Z + 2) 2 - (\ mi \ - \m 2 \) 

(Z + l)(i + 2) 3 (i + 3) 
[c 2 + .g 2 (/ 2 - 4)] 2 [Z 2 - (K| + |m 2 |) 2 ][? 2 - (Kl - |m 2 |) 2 ] 
Z 3 (Z 2 — 1) 



(4.34) 



Using the same iterative procedure one can obtain the angular eigenvalue B to higher orders in a\ and a 2 , where 
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8.g 4 w 2 [J 4 + 4Z 3 - 2l 2 (3(m 2 1 + m|) - 4) - 4Z(3(m 2 + m|) - 2) + (5(m 2 - ml) 2 + 8(m 2 + m|) - 16)] 

(Z 2 — 4)(Z - 1)Z(Z + 3)(Z + 4) 
[c^ 4 - g 2 ^ 2 (3l 2 + 12/ + 20) - 4gH(l + 4)][c^ 2 + g 2 l(l + 4)][(Z + 2) 2 - (jmij + |m 2 |) 2 ][(Z + 2) 2 - Qmij - |m 2 |) 2 ] 

Z(Z + l)(Z + 2) 5 (Z + 3)(Z + 4) 
[-w 4 + ff 2 w 2 (3Z 2 + 8) + 4. 9 4 (Z 2 - 4)}[uj 2 +g 2 (l 2 - 4)}[l 2 - (\ mi \ + \m 2 \) 2 ][l 2 - (\m x \ - \m 2 \) 2 



- 1)(Z 2 - 4) 

This is the main result of this paper, 2 where some values are presented in Figure 1 andTable I. These are compared 
with an exact numerical procedure developed in the next section. 

Note, this power series expansion of B could be used to analyze the radial scalar perturbation equation for a general 
five-dimensional Kerr-(A)dS spacetime. As of yet only the case of g = has been considered for two rotations, e.g., 
see [6, 18] (however, see [9]). 



2 A Mathcmatica notebook with results up to 8th higher order is available at http: //www-het .phys . sci . osaka-u. ac . jp/~naylor /AIM. html. 
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V. THE CONTINUED FRACTION METHOD 



Unlike the case for D > 6 [17], the angular equation in five dimensions has four regular singular points. In 
this situation, the equation can be transformed into Heun form and then the tried and tested method of continued 
fractions can be used [25]. Note, the asymptotic iteration method could also be applied [17, 26] to solve for the angular 
eigenvalue. 

To get a continued fraction we first transform the angular equation (3.6) according to 

fl(C) = C^(C-i)^(C-*o)%(C) (5.1) 

where C = (cos(20) + l)/2 and 

1 - g 2 a 2 2 

ZQ = 72772 727 ■ ( 5 - 2 ) 



Then the angular mode equation can be written in Heun form [25]: 

d 2 ( '7 S e \ d a(3( — q 



y(0 = o , (5.3) 



with the constraint 



dC 2 VC C-i (-z JdC C(C-l)(C-«b)J" 

a + (3+l=^ + 5 + e, (5.4) 



where 



a = ^(l m il + l m 2| + w/ff), P = ^(|m-i| + |m 2 | + w/g) + 2 7 = |m 2 | + l, 5=|mi| + l, e = w/g+l, (5.5) 

and 

1 oj 2 - Bg 2 m\ 1, . , 1 - g 2 a\ ,0,2, w 

1 = T~Z — 2 2 r + 7v r " 2 + ^ 9)( m 2 + ^/.9 + 2) - - 2 2\{u g - (mi + m 2 mi + m 2 + 2)), 

4g 4 af — a 2 4 4 \g z (a\ — a 2 ) 

(5.6) 

A method identical to this was used by Kodama et al. [15] for the case of a singly rotating black hole in AdS. A 
key point of the Heun form is that is satisfies a three term recurrence relation: 

a ci + p co = (5.7) 
a p c p+ i + P p c p + 7pC p _i =0, (p=l,2, ...), (5.8) 



where 



(p+l)(p + r-a + l){p + r- f3 + + 5) 
(2p + r + 2)(2p + r + 1) 



«p = /o„ 1 „ 1 n\/0^ I „ I IN ' ( 5 - 9 ) 



ep(p + r)( 7 -(5) + [p(p + r) + a/?][2p(p + r)+7(r-l)] 

ft = 7^— , 1W0 , 7\ z p(p + r) -q , 5.10 

{2p + r + l){2p + r — 1) 

(p + a-l)(p + j 8-l)(p + 7-l)(p + r-l) 

7p " (2p + r-2)(2p + r-l) ' 



with 



r = |mi| + |ma| + l. (5.12) 
The eigenvalue B can then be found (for a given ui) by solving a continued fraction of the form [27, 28]: 

a 7i «i72 "273 



ft- ft- ft- 



.. = . (5.13) 



We should finally mention that the limit g = cannot be taken with the above method, because of the divergence 
in the exponent of equation (5.1). However, if the g — > limit is taken from the start it can be considered as a 
separate case, because it removes one of the regular singular points (from four to three) and reduces to a normal 
continued fraction solution. Incidentally, the asymptotic iteration method can be used even without the additional 
factor (C — z§yl 2g in the scaling and one can find the angular eigenvalues at g = (see [26] for a 2 =0). 
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TABLE I: Comparison between the small ai,a2 expansion and the exact result for I = 0, mi = 1,7712 = for given values of 
e = ai/a2, w' = a\ui = 1 and g' = dig — 0.01. In the examples given below we found agreement between the numerical method 
and 6th & 8th order accurate to 6 s.f. 



(e, ai) 


(0.1, 0.5) 


(0.1, 1.0) 


(0.5, 0.5) 


(0.5, 1.0) 


(0.9, 0.5) 


(0.9, 1.0) 


Numerical 


3.16715 


3.66516 


3.18727 


3.74711 


3.23408 


3.93621 


2nd order 


3.16743 


3.66973 


3.18743 


3.74972 


3.23409 


3.93637 


4th order 


3.16715 


3.66521 


3.18727 


3.74711 


3.23408 


3.93621 



VI. RESULTS & DISCUSSION 



In summary we have separated the Klein-Gordon wave equation on a five-dimensional Kerr-(A)dS background with 
two rotation parameters and derived the radial and angular equations. The emphasis was on the scalar spheroidal 
harmonics, which are of relevance to charged and non-charged Kerr-(A)dS solutions [9]. We developed perturbation 
theory and have presented the result for the angular eigenvalue, B, up to 6th order (see footnote 2). This now fills 
the gap, since a perturbative expansion for the angular eigenvalues in D > 6 has been studied in [17]. 

In Figure 1 we show some plots for for varying oj',g' for fixed values of mi and m,2- As expected we find improving 
agreement as the order of expansion increases. Note that in the left panel the 2nd order answer is only linear because 
for I = mi = 7?i2 = at 2nd order there is no g dependence, cf. equation (4.36). Table I also shows data for the 
I = 0, mi = 1,7772 = mode, compared to the continued fraction method accurate to a precision of 6 s.f., which 
further confirms the improvement in increasing order. 

It may be worth mentioning that we might also be able to form a perturbative expansion using invert continued 
fractions; this was successfully applied to simply rotating cases (0,2 = 0) in [26, 28], where the natural choice appears 
to be in powers of oj' = a\ui and a,\ — a 2 g 2 , cf. [26]. However, in the case of two rotation parameters there are 
clearly two choices: one is an expansion in powers of oj' = aiu and ai — a 2 g 2 , a2 = a^g 2 ', the other (used in this 
paper) is just an expansion in powers of a\ and 02- Preliminary work on an expansion of the inverted fraction for 
two rotations seems difficult to apply symbolically. However, it is possible to verify 3 that a perturbative expansion in 
terms of Lo',a±, U2 (albeit more cumbersome) does agree with the inverted fraction second order answer [26, 28] for a 
simply rotating black hole (a 2 = 0), when a 2 ,Q!2 — > 0. We leave these interesting issues for future work. 

In conclusion, given an analytic expression for the angular eigenvalue, this can now be used to find QNMs or 
absorption probabilities (greybody factors) for a massless scalar field separated on a Kcrr-(A)dS background. As we 
mentioned in Section I this also applies to charged Kerr-(A)dS solutions (e.g., see [9]) and although we have not shown 
it we expect that this is also true for the Kerr-(A)dS-NUT case. Interesting future work might be to look at QNMs 



3 H. T. Cho, unpublished notes. 
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for general five-dimensional solutions, particularly given their relevance to gauged supcrgravity models. 
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